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Abstract 

Cellulose is one of the most versatile substances in the world. Its immense variety of applications 
was in recent years complemented by nanotechnological applications such as cellulose nanoparticle 
dressed surfaces for filtration purposes or cellulose matrices for microelectronics. The fabrication of 
such complex materials asks for thorough understanding of the surface structure and its interactions 
with adsorbates. In this study we investigate several surface model systems of nanotechnological 
interest, which are obtained by reorganization of the cellulose- vacuum or cellulose- water interfaces 
of slabs of crystalline cellulose. To do this, we equilibrated first bulk supercells of different cellulose 
allomorphs, which were constructed from crystallographic data, and then optimized the interface 
structures. From the bulk and surface systems we calculated structural properties such as unit 
cell parameters, dihedral conformation distributions, density profiles and hydrogen bonding. The 
results suggest that no overall geometrical restructuring occurs at the interface. However, the 
hydrogen bond network is strongly reconstructed, as is inferred from the dihedral conformations 
and hydrogen bond occurrences, although only within the first few layers. This holds for low index 
close packed structures as well as for high index loosely packed surfaces. Replacing the vacuum 
by ambient pressure water molecules we find less rearrangements of the cellulose surface, because 
the water allows formation of hydrogen bonds similar to those in the bulk phase. The water near 
the cellulose surface shows, however, strong structural changes. We observe reduced mobility of 
the water molecules, which corresponds to a cooling of water by about 30 degrees, in a slab that 
is about 10 A thick. Although structuring and adsorption is observed on all surfaces, no actual 
penetration of water into the cellulose structure could be observed. This suggests that pure water 
is not sufficient to produce cellulose swelling at mesoscopic timescales. This work lays the basis for 
current quantum chemical investigations on specific interaction terms within cellulose. 
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I. INTRODUCTION 

Cellulose is one of the most abundant materials on earth Exploitation of the vast 
amount of this renewable substance dates back until ancient times. The recent uprising 
of nanotechnology has opened a new range of applications for cellulose and its derivatives. 
It is used for cellulose nanocomposites, such as flexible, transparent matrices for organic 
light emitting diodes (OLED) jsj, or stable cellulose nanoparticles jsj]. To obtain new smart 
materials the surface properties of crystalline cellulose must be properly modified. The 
importance of this branch of carbohydrate research has been underlined by approval of a re- 
search proposal within the European Union 7th framework program, called SURFUNCELL, 
wherein 15 academic and industrial partners try to formulate nanostructured cellulose coat- 
ings with a variety of proposed applications, e.g. self cleaning surfaces or waste water 
treatment. 

A systematic formulation of technologically interesting materials such as the one planned 
in the mentioned project needs deep understanding of the underlying molecular structure of 
and reactivity at surfaces of the solid phases. This is, however, complicate because cellulose 
is a polymorph with several crystalline forms such as la, 1/3, II, III, but also amorphous 
forms which are partly coexistent. Cellulose is built from (3(1 — > 4) glucan molecules 
which are linear chains of D-glucose monomers linked by glycosidic bonds. In the two 
allomorphs la and 1(5 of native cellulose all chains are parallel aligned having the same 
growth direction, the allomorphs differ in their crystal structures: la has a triclinic unit 
cell of PI symmetry with a single cellobiosyl unit as basis, the crystal consists of parallel 
layers of one type, 1/3 has a monoclinic unit cell of P2\ symmetry with two cellobiosyl 
units as basis, termed origin and center, accordingly the crystal is made of two different 
alternating sheets. Regenerated cellulose II has a monoclinic unit cell of P2\ symmetry 
with two cellobiosyl units. |4j] Cellulose III is obtained when either cellulose I or cellulose 
II is treated with liquid ammonia, the resulting allomorphs are consequently denoted as 
IIIx or III 2 , both allomorphs are believed to have P2\ unit cells. Jsj] This transformation is 
reversible. Heating of cellulose III leads to cellulose IVi and IV2, respectively, which can 
reversibly be transformed into cellulose I and II, respectively, both allomorphs have a P\ 
unit cell. 

The fact that native cellulose is composed of crystalline and amorphous domains makes 
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the structural investigation very complicated. Nonetheless, many different X-ray structures 



of cellulose were published over the years j^. I Ol 1 1 Ol] . These reports gave a very fundamental 
understanding of the molecular structure, but many structural aspects are inaccessible by 
X-ray diffraction, such as the hydrogen bond (HB) network or the structure of cellulose 
surfaces or interfaces between different phases. 



Molecular modeling and especially molecular dynamics (MD) represent reliable tools 
to investigate structural and dynamical properties based on preliminary experimental re- 
sults. Considering polysaccharides and in particular cellulose many different modeling 



approaches are presented in literature, reaching from studies of bulk properties llMl5|. 



thermal res pon se 



16], relative stability 17, 18] and deformation Jl9 | , t o studies of solvated 



crystallites |2CH22j . surfaces properties or interaction with water [23N26I] and adsorbates|27|. 



Films of regenerated cellulose are semicrystalline or amorphous, that means that highly 
ordered crystalline domains of cellulose II coexist with disordered domains within their ele- 
mentary supramolecular units. Crystalline and amorphous cellulose differ in many physical 
and chemical aspects, amorphous cellulose is, for example, accessible to water, in contrast 
to crystalline cellulose. 28j Since most chemistry with cellulose is done in aqueous solutions, 
understanding of surface processes like adsorption needs a detailed knowledge of properties 



of the cellulose-water interface, however, 
structure of cellulose films are available. 



ittle data on the crystalline ordering and crystal 
29[ Modeling of the heterogenous surface struc- 
ture of cellulose films is, therefore, rather difficult, it is more difficult for the amorphous 
than for the crystalline domains. When the crystalline domains are treated as crystallites 
the surfaces are best described by crystal planes; modeling of amorphous domains that are 
accessible to water can in principle be done with very large supercells and very irregular 
surfaces, but there are no data available, which can be used for modeling these surfaces. 
It is assumed that initial solubilization processes with water are only possible when many 
polar groups are exposed at the surfaces and when the surfaces area is maximal. In order to 
investigate surface reorganization of cellulose we choose surfaces generated from large Miller 
index crystal planes and compare them with low index surfaces. We specifically address the 
reconstruction of the cellulose/vacuum and cellulose/ water interfaces and the influence of 
the surface morphology on the properties of the water component. 



II. METHODS 



Throughout this study, molecular dynamics (MD) and optimization runs were performed 
with the DL_POLY 2.20 package 



atom force field GLYCAM06 
and a non-rigid TIP3P model 
used Amber Tools 11 



30| using periodic boundary conditions. We used the all 
ycam06f.dat) [3 1] to describe the cellulose crystal systems 
321 ] to describe the water molecules. For the system setup we 



(G 



331 ] and ChemShell 3.3.1 [341 ] . therewith producing DL_POLY readable 



brce field files. Manipulations on the molecules were done using the model builder Aten 



35] , and visualization with VMD 1.8.6 [36J. Calculations were performed on an AMD64 
Linux 24 node cluster with double precision. 

All surface models are derived from the respective optimized bulk models, therefore, 
we first optimized the bulk properties of the investigated allomorphs. Each bulk model is 
represented by supercells which are repeats of the respective unit cells. The la supercell has 
5x4x5 repeats of the triclinic unit cell and the 1/3 supercell has 2x5x5 repeats of the 
monoclinic unit cell along the corresponding a, b and c directions. The cellulose II (010) 
surface model is based on a bulk model which has 5x2x5 repeats of the monoclinic unit 
cell, termed Ha. Unit cell parameters and coordinates of the cellobiose units were taken 
from recent X-ray diffraction studies j^- 10]. All three systems contained 4200 atoms and 
their volume was about 33 nm 3 . The cell parameters are the respective X-ray values in 
Table 1. For the cellulose II (120) surface model we constructed a second bulk supercell lib 
with 2x2x4 repeats of a unit cell with 8 cellobiose units ( a = 29.27 A, b = 8.93 A, 
10.36 k,a = P = 90°, 7 = 80.15°) 5376 atoms in a volume of 44 nm 3 . (cf. Fig. [1]) 

These systems were equilibrated and optimized in sequences of molecular dynamics and 



geometry optimizations, similar to a simulation protocol of [13j: 20 ps NVT MD at 600 K, 
geometry optimization, 20 ps NPT MD at 450 K, geometry optimization, 200 ps NPT MD 
at 300 K, geometry optimization. For the bulk models of cellulose II this sequence was 
preceded by a 40 ps NaT MD at 300 K and geometry optimization. 

The surface models are represented by slabs which consist of the corresponding optimized 
bulk supercells augmented by a vacuum slab of the same size. To model the la (100), 1/3 
(100), II (010) and the II (120) surfaces we thus doubled the size of our bulk supercells in one 
direction which was the b, a, b and b direction, respectively. All modeled systems consist of 
4 sheets of cellulose chains parallel to the respective surfaces, with the lowest cellulose sheet 



(200) 



FIG. 1. Schematic representation of a cellulose I/? surface. The small parallelogram at the left 
bottom depicts the unit cell of the system. The sheets are separated by dashed lines and defined 
as slabs of cellulose chains orthogonal to the envisaged surface (represented by the black line at 
the top). In this case the sheets defined by the elementary cell also act as layers in the slab. 

being frozen to simulate a bulk boundary. To find out, whether this choice of slab is a too 
severe size limitation we performed some test calculations on cellulose 1/3 and Ha with up 
to 8 sheets but did not find significant differences with respect to the measured observables, 
therefore the system size was restricted to 4 sheets throughout this study. 

The surf ace/ vacuum systems were equilibrated by the sequence: 20 ps NVT MD at 600 K, 
geometry optimization, 200 ps NPT MD at 300 K, geometry optimization. 

For the surface/water systems we started from the respective optimized surf ace/ vacuum 
systems where the vacuum slab was replaced by a 22 A thick layer of water (density p = 



1 g cm 3 ). The layer of water was built using the tool Packmol[37]. The systems were treated 
as follows: geometry optimization, 500 ps NVT MD at 300 K, geometry optimization. 

All molecular dynamics runs were performed using a 1 fs time step, Ewald summation 
and Nose-Hoover thermostats and barostats with relaxation times Tt — 0.1 ps and r p — 1 ps. 
Observed quantities for bulk and surface systems were averaged over the last 100 ps of the 
corresponding 300 K MD run. 

The self-diffusion coefficients for water were calculated using the Einstein relation. The 
mean-squared displacement data were averaged over the last 200 ps of the trajectory. The 
short averaging time does not allow an error estimate, but averaging over a significantly 
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FIG. 2. GLYCAM06 naming convention used in this work. 



longer run did not show any qualitative changes to the results. 
In this work we use the GLYCAM06 

ED. 



311 ] naming convention for glucose residues (Fig. 



III. RESULTS AND DISCUSSION 



Bulk systems 



Starting from the experimental structures obtained from X-ray diffraction, we equili- 
brated the bulk systems according to the scheme described in section [TTJ Only few of the 
relaxed cell parameters listed in Table [J differ considerably from the experimental cell pa- 
rameters. We find an increase of the cell volume by about 5% with respect to the X-ray 
structures, which is mainly caused by the much larger c parameter, whereas the a and b 
parameters as well as the cell angles are very close to the experimental starting values. The 
reason for the larger c parameters is probably the underestimation of the puckering of the 
glucose rings by the force field resulting in natter pyranose rings and, therefore, longer glu- 
can chain. This seems to be a general problem of force fields as comparison with the results 
of other modeling studies for cellulose 1/3 shows (see Tabled]). 



A very recent study by Zhang et al. 38( investigated, using MD simulation in the NPT 
ensemble and the GLYCAM06 force field, the structural response of crystalline 1/3 to heating 
from room temperature to 600 K. Experimentally, the change from a low temperature struc- 
ture to a high temperature structure is found to be reversible, but in their MD simulation 



TABLE I. Unit cell parameter of NaT equilibrated systems (this study), and values from the 
literature. Experimental studies are marked with Exp. 



Crystal phase 


Study 


a [A] 


b[A] 


c[A] 


a [deg] 


P [deg] 


7 [deg] 


V[A 3 ] 


Cellulose la 


This Study 6.89(3) 5.83(3) 10.75(1) 120.4(4) 110.6(4) 


81.7(3) 


348.3 




[10] (Exp) 


6.717 


5.962 


10.400 


lis ns 
llo.Uo 


11/18 
114. o 


80.37 


333.3 


Cellulose 1(3 


This Study 7.90(3) 


8.40(1) 10.75(1) 


8^ 1 (1} 


on i y(A\ 


102.2(5) 694.6 




M (Exp) 


7.784 


8.201 


10.380 


on 
yu 


on 
yu 


96.5 


658.3 




[13] 


8.378 


8.168 


10.523 


80 07 
oa .y I 


80 07 


90.92 


720.0 




[20] 


8.47 


8.11 


10.51 


90 


90 


90 


721.9 




[22], 298 K 


7.63 


8.23 


10.80 


89.99 


89.99 


97.17 


672.9 




[22], 500 K 


8.11 


8.28 


10.78 


90.00 


89.96 


98.33 


716.2 


Cellulose Ha 


This Study 


8.02 


9.32 


10.77 


90.2 


91.6 


121.5 


686.1 




[4J (Exp) 


8.01 


9.04 


10.36 


90 


90 


117.1 


667.8 




[8] (Exp) 


8.10 


9.05 


10.31 


90 


90 


117.1 


672.8 


Cellulose lib 


This Study 


8.30 


8.77 


10.73 


89.7 


91.2 


117.1 


694.1 



the phase transition was irreversible. The parameters of the unit cell for both structures in 
Table [I] show mainly an increase in the a parameter and a slight increase of 7, giving an 
increase of the volume of the unit cell. 

Chen et al. 



391 ] found the same irreversible change from the low to the high temperature 
structure in a simulation with the united atoms force field GROMOS 53a6 and attributed 
this fault to the parameters of the torsion potential for the exocyclic hydroxymethyl group. 
With properly changed potential parameters the phase transition became reversible. 

There is no full agreement between our cell parameters for 1(3 and those by Zhang et 
al, but this is not surprising, because Zhang et al. systematically heated the system from 
room temperature to 500 K, whereas in the mixture of MD runs and geometry optimization 
described above for the equilibration of the bulk systems the highest temperature was 600 K. 
Moreover, the super cell we used had different size and shape. Our a parameter is in-between 
the corresponding low and high temperature values of Zhang et al., the c parameters agree 
very well as do the (3 values. Our b value is closer to the high temperature value, our a 
value is by 5 degrees smaller than both values by Zhang et al., and our 7 value is 4 and 
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FIG. 3. Different conformations of C6 group, x is denned as C4-C5-C6-06 dihedral angle. 



5 degrees larger than the high and low temperature values by Zhang et al, the volume of 
our elementary cell is in-between the low and high temperature values. We assume that the 
equilibrated bulk cell obtained with our scheme is neither a full low temperature nor a full 
high temperature structure. 

The phase transition is accompanied by a change in the conformation of the exocyclic 
hydroxymethyl group, which can be found in three staggered conformations termed gg, gt 
and tg where the first g and t mean gauche and trans position, respectively, of the C6-06 
group with respect to C5-05, whereas the second g and t mean gauche and trans position of 
the C6-06 group with respect to C5-C4. The conformations can be characterized by, e. g., 
the values of the C4-C5-C6-06 dihedral angle x (Figure E]). In the tg and gt conformations 
the position of the C6-06 bond is equatorial with respect to the glucose six ring, in the 
gg conformation the C6-06 bond is axial. In the tg conformation can accordingly the 
06H group participate in intra- and inter-chain hydrogen bonds, in the gt conformation it 
participates only in interchain hydrogen bonds and in the gg conformation the 06H group 
can participate only in inter-sheet hydrogen bonds. Zhang et al. found that upon heating tg 
conformations in the origin sheets are dominantly changed into gt conformations, but into 



gg conformations in the center sheets. [381] Chen et al. [39j showed that above 500 K the tg 
conformations are completely replaced by gt and gg conformations. 

Figure H] shows the distribution of the hydroxymethyl conformations in the 4 cellulose 
systems investigated. In 1(3 the ration of tg : gt : gg is 38 : 31 : 31, so our \f3 system 
is neither the low temperature system with 100 percent tg, nor is it the high temperature 
system with no tg conformations at all. The equal population of gt and gg conformations 
is in accord with the findings of Zhang et al. and Chen et al. Considering the influence 



a 




gg gt tg 



FIG. 4. Frequency of the hydroxymethyl conformations. 

of the force field on the population of the hydroxymethyl conformations, one can assume 
that the high gt population in la is also an artefact of the force field, the extremely low gg 
population is in agreement with experimental findings. 

For the investigation of cellulose II we used two different supercells, the distribution of 
the hydroxymethyl conformations differs considerably for the gt and tg conformations: in 
both systems tg dominates but in Ha the percentage is below 50 percent and in lib it is 
about 66 percent. In both systems the gg population is about 20 percent. We assume that 
the difference is caused by the equilibration scheme and the different size and shape of the 
supercells. 

Hydrogen bonds in cellulose bulk 

A hydrogen bond (HB) is a system where a hydrogen atom is covalently bonded to a 
donor oxygen OD-H bond, which weakly interacts with the lone pair electrons of an acceptor 
oxygen atom OA, that is not directly bonded to the group. Such a HB is represented by 
ODH...OA. Since the OH bond lengths in different HBs are rather constant, the non-bonded 
O...H distance and the 0H...0 angle are sufficient to characterize a HB. 

Two types of intra-chain HBs between glucose rings can be found in a glucan molecule; 
the first one is 03-H...05 with the hydroxyl group 03H in one glucose ring and the acceptor 
atom 05 in the next glucose ring; the second HB is 02-H...06 with the hydroxyl group 

10 



TABLE II. Distribution of hydrogen bonds in la bulk together with HB parameters. Only hydrogen 
bonds with a percentage higher than 0.5 are shown. Distances are given in Angstrom, angles in 
degree. 





HB 


type a rel. 


occ. d(0 ■ 


■HO) Z(OHO) 


05- 


• • H03 


A 


32 


1.86 


155 


06- 


••H02 


E 


32 


1.88 


153 


06- 


■H03 


E 


1 


2.06 


159 


03- 


••H06 


E 


16 


1.91 


157 


02- 


••H06 


A 


16 


1.93 


159 


02- 


••H03 


A 


1 


2.09 


146 


04- 


•H06 


S 


1 


2.01 


153 



a) A means intra- and E means inter-chain, and S means inter-sheet hydrogen bonds. 

02-H in one ring and the acceptor atom 06 in the second ring. This HB is only possible 
when the exocyclic hydroxymethyl group has tg conformation. Thus, the covalent glycosidic 
bond and the two HBs build a strong link between two glucose rings and make the chain 
rather stiff. Interchain HBs connect glucan molecules to sheets, which interchain HBs are 
formed depends on the crystal structure and on the temperature. 

Table [TT] shows the relative occurrence of the different HBs in the equilibrated la bulk. 

Most important are the intra-chain HB 03H...05 and the inter-chain HB 02H...06; with 
32 percent each, they represent nearly two third of the hydrogen bond network in the bulk, 
to the remaining third contribute the intra-chain HB 06H...02 HB and the 06H...03 inter- 
chain HB with about 16 percent each. That means, about half of the hydroxymethyl groups 
are found in the tg and the other half in the gt conformation. The rest of the HBs are weak 
intra-chain and inter-sheet HBs. 

Crystalline cellulose 1/3 is made of two alternating sheets (Figure EJ) corresponding to the 



two different cellobiosyl units in the unit cell, called center and origin units. In one shee ; 
chains are nearly coplanar, in the other sheet the chains are slanted, Heiner et al. 23, 



the 



40| 



called these sheets "odd" and "even", we prefer "flat" and "rippled". According to authors 



like Nishiyama et al. IQj the origin units create the flat sheets others like Zhang et al. [38 1 
make the flat sheet with the center units. Nishiyama et al. js] showed for the first time the 
existence of different HB networks in the two sheets, Figure [5] visualizes the different HB 
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pattern. 




FIG. 5. a,b base plane projection of the cellulose 1/3 (100) slab. 

The analysis of the HB network in equilibrated 1(3 bulk shows that in the rippled sheets 
there are about 20 percent more HBs than in the flat sheets, which corroborates the existence 
of two different HB networks. A detailed comparison of the two networks reveals equal 
numbers of intra-chain HBs 03H..05 and inter-chain HBs 02H..06 in both sheet types, 
these two HBs form the equal backbone of the HB networks. The difference between the 
two sheet types comes from the different numbers of 06H..03 inter-chain HBs, which is in 
the rippled sheets about 2 to 3 times higher than in the flat sheets. That means, the stiffer 
chains in the flat sheets are connected mainly by 02H..06 inter-chain HBs, the less stiff 
chains in the rippled sheets are connected by the same number of inter-chain HBs 02H..06 
and by additional 06H..03 inter-chain HBs. 

Compared with la, the number of tg conformations of the hydroxymethyl group decreases 
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TABLE III. Distribution of HBs in 1/3 bulk together with HB parameters. Only bonding atoms 
with a percentage higher than 0.5 are shown. 



HB 


rel. 


occ. 


d(O---HO) Z(OHO) 








bulk 




03H- ■ 


05 


34 


1.83 


157 


02H-- 


06 


34 


1.91 


158 


03H- • 


06 





2.07 


148 


06H- • 


03 


14 


1.99 


154 


06H-- 


02 


12 


2.00 


154 


06H-- 


04 


6 


1.96 


157 








flat, slippts 




03H- • 


05 


42 


1.83 


156 


02H-- 


06 


40 


1.96 


156 


03H-- 


06 





2.01 


150 


06H- • 


03 


11 


2.05 


157 


06H-- 


02 


7 


2.01 


155 


06H- • 


04 





2.09 


138 






rippled sheets 




03H- • 


05 


35 


1.83 


157 


02H-- 


06 


34 


1.91 


161 


03H-- 


06 


1 


2.04 


149 


06H-- 


03 


23 


1.96 


152 


06H-- 


02 


8 


2.14 


155 


06H-- 


04 





2.03 


137 



in both sheet types similarly strong, but the increase of gt and gg conformations in the 
two sheet types is very different: in flat sheets there are more gg conformations than gt 
conformations, in rippled sheets there are nearly no gg conformations. Consequently, most 
hydroxymethyl groups in rippled sheets are involved in inter-chain HBs, whereas in fla t 
sheets they are also involved in inter-sheet HBs. This was also found by Zhang et al. 
Relative frequencies and geometry parameters of the HBs are given in Table IIHI 
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This result diverges considerably from the C13 CP/MAS results that the hydroxymethyl 



411 ] Experimentally, two different hydrogen 



group in la and 1/3 has only the tg conformation, 
bond networks exist in 1/3, but in both networks the hydroxymethyl group has only the tg 
conformation. 



In analogy with cellulose 1/3, monoclinic cellulose II crystals are made of alternating 
sheets, but the growth directions of the glucan molecules are reversed in alternating sheets. 
This allows HBs which do not occur in 1/3 and, thus, gives rise to different topologies of 
the HB network. Inspection of the a,b base plane projection of the cellulose II supercell in 
figure [6] shows clearly the existence of two sheet types but a distinction between flat and 
rippled seems to be rather artificial. In the a,b base plane projection of the origin sheets, 
which are made of the origin units in the unit cell, one can observe a cyclic structure of the 
inter-chain HBs (sheets 1 and 3 in Figure [6]), we find mainly 06H..02 and little 06H..03 
HBs, the structure is highly regular as nearly all 06H groups are involved in these HBs; the 
alternating position of the glucose rings that are connected by the inter-chain HBs causes the 
regular HB network. That means also that in origin sheets there are no 06H..02 intra-chain 
HBs; the intra-chain HB 03H...05 has the same frequency as in 1/3. 



For the center sheets, which are made of the center units, no cyclic structure of the 
inter-chain HBs can be seen in the a,b base plane projection. We observe a low regularity 
in the HB network, which is mainly made from 02H..06 and 06H..03 inter-chain HBs, 
02H...06 HBs occur nearly exclusively in the center sheet. Not all 06H..02 intra-chain 
HBs are broken and replaced by inter-chain HBs, but the number of 03H...05 intra-chain 
HBs is in the center sheets only half of that in the origin sheets, the 03H groups are mainly 
targets of 06H groups. In center sheets one can also find 02H..06 intra-chain HBs. 



Inter-sheet interaction mainly happens via 02H..02 HBs with the 02H group from the 
origin sheet, and 06H..02 HBs, with the 06H group from the center sheet, but there are 
also 03H..06 and 06H..06 inter-sheet HBs. 



Note, the relative occurrences for the bulk in Table HVl are not the sum of the occurrences 
for the layers, the latter do not contain any inter-sheet HBs. 
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TABLE IV. Distribution of HBs in Ha bulk together with HB parameters. Only bonding atoms 
with a percentage higher than 0.5 are shown. 



HB 


rel. 


occ. 


d(O---HO) Z(OHO) 


bulk 


03H- • • 




23 


1 Q4 


153 6 


02H- • • 


06 


17 


1.82 


161.3 


06H- • • 


03 


17 


1.85 


156.2 


06H- • • 


02 


14 


1.83 


158.8 


03H- • • 


06 


10 


2.03 


153.0 


06H- • • 


06 








06H- • • 


04 








02H- • • 


02 


17 


1.85 


2.830 












03H- • • 


05 


47 


1.93 


155.4 


02H- • • 


06 


1 


1.98 


154.9 


06H- • • 


03 


12 


1.94 


160.5 


06H- • • 


02 


39 


1.82 


158.9 


03H- • • 


06 


1 


1.82 


157.6 


06H- • • 


06 


1 


1.98 


154.9 


06H- • • 


04 














center sheet 




03H- • ■ 


05 


21 


1.96 


149.4 


02H- • • 


06 


47 


1.81 


161.0 


06H- • • 


03 


28 


1.82 


155.0 


06H- • • 


02 


1 


2.12 


150.0 


03H- • • 


06 


2 


2.02 


150.7 


06H- • • 


06 








06H- • • 


04 





1.97 


139.3 
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Thermodynamic stability of cellulose bulk 

By averaging over configurational energies of the equilibrated systems we get an estimate 



of the relative thermodynamic stabilities of the different allomorphs. From experiment 171, 



421 ] it is known that heating of cellulose la will convert it into 1/3 and finally into II, but not 
the other way around. Experimentally, cellulose II is thermodynamically most stable and 
la is least stable, our simulations yield the following ordering of the energies per cellobiose 
unit: 215.0 kJ mol" 1 (1/3), 213.8 kJ mol" 1 (la), 211.6 kJ mol" 1 (Ila) and 210. 7kJ mol^ 1 
(lib). 

Our results corroborate the finding that cellulose II is more stable than cellulose I, the 
d(0 ■ ■ ■ HO) distances are in cellulose II shorter than in 1/3, indicating that the glucan 
molecules are closer together and this will increase the stability because the HB and the 
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FIG. 7. Schematic structure of cellulose la (a), 1/3 (b) and II (c) nanocrystallites, exposing different 
crystal planes. 

dispersion interaction will be stronger. According to our quantum chemical calculations 
dispersion interaction between the chains will be more important than HBs. The expected 
higher stability of cellulose 1/3 with respect to cellulose la was, however, not found. We 
attribute this discrepancy to the fact that our equilibrated 1(3 system is indeed a less stable 
high temperature system. We are sure that the difference in the stability of cellulose Ha and 
lib is caused by the equilibration scheme. Since we were mainly interested in the generation 
of surfaces of different roughness for the study of the cellulose-water-interface we did not 
investigate the origin of the difference between the two supercells. 



Cellulose/ vacuum and cellulose/water interfaces 

The aim of this study is to investigate how the properties of the cellulose/water interface 
depend on the structure of the cellulose surface. Macroscopic cellulose samples, such as fibres 
or films, are often regarded as being composed of nanocrystallites embedded in amorphous 
cellulose; in a first approximation, the surface of such a sample can then be seen as a tiling 
made of many different crystal faces of irregular shape, which are described by the Miller 
indices of the corresponding crystal planes (see Figure [7]) . As models for surfaces of the 
amorphous parts we choose crystal planes with large Miller indices showing large surface 
roughness. The properties of the macroscopic surface is then an average of the properties of 
all exposed crystal faces. 

It is not possible to investigate the surface properties for all exposed crystal planes, but 
only for few representatives which differ most strongly in certain properties like surface 
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roughness, polarity etc. and in the interaction with adsorbates. This approach has also 
been chosen by e. g. Mazeau and Vergelati . In the la (100) system the bulk sheets are 
not parallel to the surface, in the layers parallel to the surface the glucose rings are slanted; 
in the 1/3 (100) the layers are parallel to the surface and at least in the flat sheet also the 
glucose rings are parallel; in II (010) the layers are made of parallel sheets but the glucose 
rings are slanted; in II (120), finally, the sheets are not parallel to the surface, the layers are 
made from glucan molecules that stem alternatively from origin and center sheets and the 
overall surface shows a zigzag behavior. 

The aforementioned molecular dynamics scheme reveals very different degrees of recon- 
struction for the different surfaces. The la (100) surface, where the rings of the cellobiose 
molecules are slanted with respect to the (100) plane, shows high irregularities after equili- 
bration (shown in figure |8b); inter-chain hydrogen bonds become inter-sheet hydrogen bonds 
and this causes a reduction of the number of hydroxy groups pointing into the vacuum and a 
strengthening of the bonding between the glucan molecules in the top layer. Not surprising 
is that the flat (100)/(010) surfaces of cellulose 1/3 and II show little reconstruction, this 
could be different for the more rippled 1/3 (200) surface, but we did not investigate it. In 
spite of the grooved structure of the cellulose II (120) surface we did not find a considerable 
reconstruction due to vacuum exposure. 

To quantify the different roughnesses of the surfaces and thereby the degree of surface 



reconstruction we adapted a measure which is very common in material science 42fl . The R a 
value of a surface is defined as follows: 



R, 



1 n 

a = - yj I Ajfcl Ay { = yi-y (1) 



n . 



Whereas y, is the absolute height of a surface atom and y defines the average surface height. 
As the average surface height we took the tangential plane to the closed van der Waals 
spheres of the uppermost layer. 

In Table |V] we show the roughnesses of the four reconstructed surfaces. Roughness is 
a property of the structure of surfaces, the values from Table [V] support the classification 
of 1/3 (100) as smooth and of II (120) as rough, that II (010) is much smoother than la 
(100) is not that obvious. Interestingly, the difference in roughness between a bulk layer, 
the surface exposed to vacuum or to a water environment is insignificant compared to the 
overall difference between the different cellulose types. In addition to the unchanged surface 
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FIG. 8. Cutout of the supercells showing a) parts of two layers of Cellulose la bulk, b) parts of 
the first two layers of a Cellulose la (100) surface exposed to vacuum and c) parts of the first two 
layers of a Cellulose la (100) surface exposed to 300 K TIP3P water. 

roughness we also were not able to find significant changes of the layer distances orthogonal 
to the surface. This suggests that at this time scale no drastic changes in the mesoscopic 
geometric structure of cellulose occurs. 

Another frequently used measure for roughness is the standard deviation of the distribu- 
tion of the absolute heights ji from the average surface height 



The two measures give different roughness values, the R a values are always larger than the 
R a values. For a large number of uniformly distributed random numbers yi the difference 




(2) 
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TABLE V. Roughness R a [ A] of cellulose surfaces exposed to vacuum. 

Surface R a [A] 
Cellulose la (100) 0.7418 
Cellulose 1(3 (100) 0.5429 
Cellulose II (010) 0.5890 
Cellulose II (120) 0.8289 



is about 16%, for normally distributed random numbers the difference is about 25%. The 
coordinates of the top layer atoms are neither uniformly nor normally distributed random 
numbers, therefore, we regard the larger roughness values R a reported by Mazeau as 
consistent with our data. 

Nonetheless, surface reconstruction at the solid/vacuum interface is accompanied by a 
significant change of the molecular hydrogen bond network, which can be described by the 
frequencies of occurrence of the three conformations of the exocyclic C6 group and by the 
number of hydrogen bonds. The data were obtained from the last 100 ps of the 300 K MD 
runs. 

In Figure [9] the dihedral angles of the exocyclic C6 group are plotted which were measured 
for the layers in la (100), but we find for all systems similar distributions with peaks that 
are easily distinguishable and can be assigned to one of the three dihedral conformations. 
The area under the peak can be interpreted as the probability of finding this conformation, 
the relative probabilities for all systems are given in Table IVII The reference data for the 
crystalline systems with equivalent layers (la and lib) are obtained by averaging over all 
layers, as for cellulose la and cellulose lib, or by averaging over the two different kinds of 
layers, as for cellulose 1(3 and cellulose Ha. The averaged data are always given in line L3, 
for the systems with two different kinds of equivalent layers, the corresponding data for LI 
and L2 are also given. 

In the cellulose la bulk system all layers parallel to the surface show the same relative 
probabilities for the dihedral conformations within small statistical fluctuations. The small 
amount of gg conformations (6 percent) and the large amount of gt (40 percent) and tg 
(54 percent) conformations once again demonstrates that in la there are nearly no inter- 
sheet HBs, but many intra- and inter-chain HBs. The alternating layers in the 1(5 system 
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FIG. 9. Representative distributions of C6 dihedral angle conformations for different layers in a 
Cellulose la (100) surface exposed to vacuum. Also shown are the integrated occurrence frequencies 
that are used for analysis. 

are identical with the flat and rippled sheets and have rather different distributions of the 
dihedral conformations with a high frequency of gg conformations in the flat sheets (layers 
1 and 3) and nearly no gg conformations in the rippled sheets (layers 2 and 4), as already 
discussed above. 

In the Ha system again the layers are identical with the sheets and we find accordingly 
different probability distributions in alternating layers. In the origin sheets (layers 1 and 
3) we find again a high number of gg conformations indicating that the C60H groups from 
these sheets form active inter-sheet HBs whereas little active inter-sheet HBs are made with 
C60H groups from center sheets (layers 2 and 4). In the lib system the layers are made 
from equal contributions of origin and center glucan molecules and, therefore, all layers 
show the same probability distribution. The difference between systems Ha and lib is again 
demonstrated in the averaged probability distributions which differ mostly in the gt and tg 
probabilities. 

Now we compare the bulk data with the probabilities for the layers of the cellu- 
lose/vacuum interfaces. For la (100) there is a higher probability to find gg and gt confor- 
mation in the layers near to the surface and a reduced probability to find tg conformations. 
The increase of gg conformers corresponds to the increase of inter-sheet HBs between the 
tilted cellulose chains at the surface. The probability distribution of layer 3 is again very 
close to the bulk distribution, we interpret this finding as showing that reconstruction of 
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the la (100) surfaces is restricted to only few layers. 

In cellulose 1/3 (100) on the other hand we find only little differences between the prob- 
ability distribution for the bulk and for the reconstructed surface, because in this system 
the layers are formed by sheets and the top layer is a flat sheet. Maybe, reconstruction of 
the (200) surface would result in larger differences at least for the top layer, but we did not 
investigate this surface. 

Although in Ha (010) the layers are identical with sheets, we find strong increase in 
the probability for gg conformations in the first two layers, the third layer shows again a 
distribution very similar to that in the bulk. In lib (120), finally, again only the two top layers 
are strongly reconstructed, the tg conformations are decreased and the gt conformations are 
increased. 

There is, thus, a very good correspondence of surface roughness and surface reconstruc- 
tion, the surface with the lowest R a value shows the least reconstruction, irrespective of 
whether the layers are sheets or not. This correspondence was also reported by Mazeau et 
al. Q for cellulose I. 

When the vacuum of the solid/vacuum interface is replaced by bulk water we find very 
different changes of the probability distributions: The changes in the two top layers of la 
are much smaller than the changes between the bulk and the solid/vacuum interface, the 
changes in layer 3 can be neglected. In 1/3, on the other hand, we see in the top layer a 
strong increase of the gg and a strong decrease for the tg conformations, in layer 2 the gt 
conformations increase in expense of the tg conformations. The changes in layer 3 are again 
small. 

In Ha replacement of the vacuum by water is accompanied by little changes in the proba- 
bility distributions, in lib a significant increase of gt in expense of tg conformations is found 
only in the top layer, the other changes are again rather small. 

Summarizing, we found that surfaces that show strong hydrogen bonding reconstruction 
at the solid/vacuum interface, change little when water replaces the vacuum and vice versa, 
and, none of the top layers resemble very strongly the corresponding bulk layer hydrogen 
bond network, that means, there is always a surface reconstruction at an interface, the kind 
and amount of reconstruction of the cellulose surface depends on the second phase at the 
interface, the cellulose allomorph and the structure of the crystal plane. Similar results 
were obtained by Newman and Davidson [45J in a 13 C NMR study of the cellulose- water 
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TABLE VI. C6 dihedral angle distributions for two cellulose interfaces and crystalline bulk. The L 
stands for the layer in the slab. In systems where the layers are equivalent, the average distribution 
for the whole slab is shown. Details are given in the main text. 



Crystal system L crystal bulk solid/vacuum solid/water 

gg : gt : tg gg : gt : tg gg : gt : tg 
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interfaces for cellulose I and cellulose II surfaces. From AFM studies of cellulose la surfaces 
Baker et a/.j^fj conclude that hydroxymethyl groups will change at the surface from tg to 
gt conformation. 

The only significant difference between results obtained from the simulations with a 4 
layer and an 8 layer supercell, respectively, was found for the distribution of the C6 dihedral 
angle in the 1/3 bulk which has an a-b-a-b-. . . layer structure. There are, however, large 
differences in the layers 1 and 3, and 2 and 4, respectively. We find 0.43 : 0.00 : 0.53 : 0.02 
for the gg conformations, 0.23 : 0.65 : 0.15 : 0.31 for the gt conformations and 0.35 : 0.37 : 
0.32 : 0.67 for the tg conformations. In the 8 layer supercell the differences between different 
flat and different rippled layers are much smaller, these data are included in table IVI( they 
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TABLE VII. Number of hydrogen bonds within layers per cellobiose unit for the bulk and surface 
systems, excluding intra-chain 03H- • • 05 hydrogen bonds. The L stands for the layer in the slab. 



Phase L crystal bulk solid/vacuum solid/water 



Cellulose la 1 


2.10 


2.28 


1.94 


2 


2.01 


2.49 


2.30 


3 


1.99 


2.30 


2.35 


Cellulose 1/3 1 


2.60 


2.92 


1.92 


2 


3.51 


3.34 


3.19 


3 


2.63 


2.48 


2.49 


Cellulose Ha 1 


2.12 


2.48 


2.11 


2 


3.01 


2.41 


2.60 


3 


2.09 


2.20 


2.22 


Cellulose lib 1 


2.26 


2.68 


2.26 


2 


2.66 


2.41 


2.23 


3 


2.36 


2.20 


2.35 



are, nevertheless larger than those found for lib, which has the same layer structure as 1/3. 

The HB networks in cellulose and at the interface are made from more than just the C60H 
groups, so we calculated also the average number of all intrachain and interchain HBs, except 
03- • • 05 HBs per cellobiose unit for each layer in the bulk and interface systems (see Table 
IVIII) . The analysis did not allow to include intersheet hydrogen bonds, but the unchanged 
layer distances suggests that the amount of interlayer hydrogen bonding is not strongly 
affected by interface reconstruction. 

Again, we find for the four systems rather different trends: in the top layer of la, that 
shows the strongest surface reconstruction, there is a moderate increase of the HBs when 
going from bulk to solid/vacuum, and a slightly stronger decrease when going to solid/water 
because the cellulose.. water HBs are not counted. In the second and third layer the number 
of HBs increases with respect to the bulk. Similar trends are found for the 1/3 system. 

In the cellulose II systems we observe in the top layer an increase of the number of HBs 
when going from bulk to solid/vacuum and a strong decrease when going from solid/vacuum 
to solid/water, there one has about as many HBs as in the crystal layer. In layer 2 of both 
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interface systems there is a decrease of the number of HBs with respect to the bulk, in layer 
3 the resemblance to the bulk is striking. 

When a surface is created by cutting a slab out of the crystalline bulk, there will be 
hydroxy groups in the top layer having no partner to form HBs between the layers (dangling 
HBs). During the reconstruction of the solid/vacuum interface dangling HBs will be partially 
converted into inter-chain or maybe also inter-sheet HBs (which where not covered in the 
analysis of table IVIip in the top layer. The additional HBs, due to vacuum exposure one can 
directly observe for Ialpha (100) in figure [HJ In the solid/water interface the former dangling 
HBs and some of the inter-chain HBs in the top layer form HBs to the water slab. This yields 
a reduced average amount of HBs within the first cellulose layer. We did not find, however, 
a significant amount of water molecules penetrating between cellulose chains and diffusion 
of water molecules into the solid bulk; the formation of hydrated hydroxy groups at the 
cellulose surface does not seem to be energetically sufficient, at least not in the simulation 
times used in this study. This leaves us to conclude that pure water at ambient conditions 
is not able to cause swelling of cellulose at microscopic timescales, but merely wetting. 

Structure and dynamics of water at the interface 

To investigate the structure of the water slabs at the interfaces we calculated density 
profiles of the systems orthogonal to the surface, which are shown in Figure [10J The black 
lines depict the density profiles of the cellulose, the red ones the density of the water. 

The cellulose profiles reflect nicely the different layering in the four systems: fine structure 
is found when the sheets are not parallel to the surface as for la and Ha, the larger the angle 
between surface and sheet is, the more structure is found in the profiles. For 1/3 and lib, 
where the sheets are parallel to the surface, we just see the peaks of the four and eight layers, 
respectively, in the systems. 

The density profiles of water have an oscillatory form indicating an ordering in the water 
slab that is not typical for bulk water, the highest density maximum is found near the 
interfaces. The deviations from the bulk density are strongest for the flat surfaces 1/3 (100) 
(40 percent) and II (010) (60 percent), and least pronounced for the non-flat surfaces la (100) 
(10 percent increase) and II (120) (20 percent increase). The oscillations decay completely 
within 10 A at the la (100) surface the water bulk is undisturbed already after about 4 to 
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FIG. 10. Density profiles of cellulose-water interfaces, the cellulose profiles (solid left) are for la 
(100), 1/3 (100), II (010) and II (120), right are the profiles (solid right). The left scale in g cm" 3 
is for the density profiles. Also shown is the diffusion coefficient of water molecules as a function 
of the surfaces distance (dashed) with the right scale in 10 -9 m 2 s _1 . 



5 A. 

When placing an undisturbed slab of water onto a cellulose slab, the water molecules 
rearrange to maximize interactions. We therefore see areas around the cellulose chains, 
where water molecules adsorb preferentially (around the terminal hydroxy groups) and areas 
where no hydrogen bonds can be formed, such as around the unpolar hydrogens bonded to 
the glucan framework. These areas show depletion of water molecules compared to the 
terminal hydroxy groups (cf. Figure [TT]). This densification and restructuring of water on 
the cellulose surfaces happens in a volume of a few A of vertical distance to the surface. Using 
a similar reasoning as 24| one can relate the degree of hydrophilicity to the density profile: 
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TABLE VIII. Diffusion coefficients in m s 2 for the first 2 A thick layer of water in contact with 
the different cellulose crystal surfaces at 300 K. 



Crystal surface 


DxlO 9 [ms~ 2 ] 


Cellulose la (100) 


1.80(1) 


Cellulose 1/3 (100) 


2.11(1) 


Cellulose II (010) 


2.30(1) 


Cellulose II (120) 


1.90(1) 



The more sharp the first peak is, corresponding to a dense layer of water adsorbed on the 
surface, the more hydrophilic is the surface. However, densification of a water slab above a 
cellulose surface can also be related to the regularity of the surface. Highly regular surfaces, 
such as cellulose 1(3 (100) and cellulose II (010), allow a large amount of water molecules 
within a slab of a certain vertical spacing above the surface. In the case of cellulose la (and 
also cellulose II (120) the tilted and slightly irregular chains distribute this effect over a 
larger vertical range above the surface, although the amount of interacting water molecules 
may be similar, if not more. The densification of a water slab above a cellulose surface is 
therefore both, a function of the hydrophilicity and the regularity (roughness) of the surface. 
We therefore refrain from extrapolating solel y fr om the density profile to the hydrophilicity 



of the surface, as has been done in literature 2J, |40], as long as we don't know in detail how 



the interaction energy at the water/surface interface depends on the surface morphology. 
The surface structure of cellulose II (120) surface with eight layers (see Figure [TO]) allows 
interaction of water molecules not only with the top layer but also with the second layer, 
this is reflected in the shifted onset of the water profile which is at the same position as the 
maximum of the density of the cellulose top layer. Cellulose II (120) should therefore have a 
larger water accessible surface area and we hoped to observe initial solvation processes. No 
such processes were observed, however, within the simulation time. This supports the finding 
that pure water and cellulose do not interact strongly enough to allow easy solubilization at 
ambient conditions. 

Another property that can be used to describe the structure of the water slab at the 
cellulose surface is the self diffusion coefficient as a function of surface distance. The non-rigid 
TIP3P water model used in this study yields for bulk water at 300 K a value of 4.0 10 9 m s~ 2 , 
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FIG. 11. a,b base plane projection of a) the cellulose la (100) and b) the cellulose 1/3 (100) slab 
in contact with water at room temperature and room pressure. Shown are the first layers of 
cellulose chains as sticks and the adsorbed water molecules in a continuous van der Waals surface 
representation. 



which agrees well with our values at 10 A, this value is about 50 percent larger than the 



experimental value of 2.4 x 10~ 9 m 2 s _1 (Refs. 47], |48J), but considerably better than the 
value that is obtained with the rigid TIP3P model (see Ref 49]), although at considerably 
higher computational cost. Disregarding the absolute values of the diffusion constant, one 
can nevertheless obtain an interesting information on the cellulose/water interface: we find 
a decrease in the diffusion constant from bulk water to water at the interface of about 
1.6 x 10 -9 m 2 s" 1 and such a decrease corresponds, according to the data by Harris et al.\^\, 
to a decrease of the water temperature of about 30 degrees. The reduction of the diffusion 
constant at the interface is caused by at least two factors: i) the reduced degrees of freedom 
of water motion, and ii) the formation of variably dense water regions ('arches') at the 
cellulose/ water interface of different stability. One could wonder if the shape of the graphs 
of the diffusion constant vs. distance from the surface reflects the differences in the water 
structures as described by the density profiles; it is however not possible to show such 
a dependence, because the smallest bin size used for the calculation of the self diffusion 
constant by integration of the velocity autocorrelation function was a slab of 2 A, which is 
larger than the strongly fluctuating interface region. In any case, the reduced water mobility 
and the formation of water structures of different stability at the cellulose surfaces, together 
with the change in hydrogen bonding patterns can be seen as a main ingredient of wet ting 



of crystalline cellulose. Our results are consistent with the results from a 2 H NMR study 
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50] 



who found three different types of water at the cellulose interface: a strongly bound, rigid 
water, which cannot crystallize, a highly mobile water, whose motion is hindered due to 
anisotropic constraints, so that it can only undergo fast 180 flips around their bisector axis, 
and finally the mobile water that can perform isotropic motions. This agrees also very well 



with the results from 1 H and 13 C solid state NMR results by Taylor et al. 5l| that moisture 
adsorbed at cellulose surfaces has different correlation times. 

As shown by Vega et al. [49], the choice of the water model crucially influences the 
properties of liquid water and water steam. It should not be much of a surprise if the water 
model also influences the dynamical properties of interfaces. Dynamical observables might 
depend strongly on how balanced the force fields used for the cellulose and water description 
are. 



IV. CONCLUSION 

Cellulose is a biological material of immense structural complexity. In this work we 
approached this material starting from idealized crystalline bulk structures. We investigated 
the structure of different bulk crystals that allows a somewhat closed view. It is often 
claimed in literature that the attractive interaction between la sheets is mainly dispersion 
interaction, whereas 1/3 and II show strong polar intersheet interactions; it should be clear 
that dispersion interaction is important for all kind of interfaces, irrespective of whether there 
are hydrogen bonds or not. Due to similar distances in our classical dynamical approach 
we always find a constant background of dispersion interaction for all crystal systems. Of 
course, one should know how much of the stabilization energy is dispersion energy and how 
much is, say, hydrogen bonding energy, but to do this quantum theoretical methods must 
be used. Work is in progress in our group to investigate the stabilization of chains in sheets, 
or between sheets in bulk, and to calculate the contribution of dispersion interaction to 
stabilize these structures. 

Knowledge of the reorganization of these bulk materials on interfaces is a basic prereq- 
uisite to understand more complex interactions that would involve ions, restructuring at 
non-standard conditions and, most interestingly, adsorption of complex nanostructured ma- 
terials. When cellulose surface slabs are allowed to reconstruct, the degree of reconstruction 
depends strongly on the allomorph and the type of crystal plane that is chosen as slab sur- 
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face. As a main result we find that within nanoseconds no significant mesoscopic structural 
changes, such as changes in layer distances, chain tilting or strong irregularities, occur. This 
might lead to the conclusion that cellulose surfaces in contact with vacuum or water show no 
restructuring at all. A detailed view on the molecular structure and the hydrogen bonding 
network reveals substantial restructuring, although only for the first few layers. 

Adsorption of water on cellulose surfaces cannot be solely explained with the formation 
of hydrogen bonds between water molecules and oxygen containing groups at the surface, as 
in the bulk, dispersion interaction must be considered as well. And to find out, whether the 
interaction between glucan molecules, that is hydrogen bonds and dispersion interaction, are 
quantitatively correct described with force fields and water models, comparative quantum 
theoretical investigations are needed, specifically addressing the hydrogen bond strength 
and amount of dispersion interaction between the chains in cellulose. Such investigations 
are under way in our group (Hoja, Grossar, Maurer and Sax, to be published). The results 
presented here point to different wettability of different cellulose surfaces, due to a different 



ratio o 
studies 



polar and unpolar moieties. This observation has already been made in recent 
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40] . We can confirm those results and augment them with the observation that 



even loosely packed high index surfaces, such as Cellulose II (120) do not allow significant 
solvation and swelling within a mesoscopic time frame. We further found hints that water 
applied to surfaces that were optimized against vacuum can partially reverse the surface 
reconstruction, even though they did not insert between cellulose sheets, that were previously 
reconstructed against vacuum. Anyhow, we did not find that swelling of cellulose can be 
caused by pure water. 
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